Magnetic response of carbon nanotubes from ab initio calculations 
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We present ab initio calculations of the magnetic susceptibility and of the ^^C chemical shift for 
carbon nanotubes, both isolated and in bundles. These calculations are performed using the recently 
proposed gauge-including projector augmented-wave approach for the calculation of magnetic re- 
sponse in periodic insulating systems. We have focused on the semiconducting zigzag nanotubes 
with diameters ranging from 0.6 to 1.6 nm. Both the susceptibility and the isotropic shift exhibit 
a dependence with the diameter (D) and the chirality of the tube (although this dependence is 
stronger for the susceptibility). The isotropic shift behaves asymptotically as a/D + 116.0, where a 
is a different constant for each family of nanotubes. For a tube diameter of around 1.2 nm, a value 
normally found in experimental samples, our results are in excellent agreement with experiments. 
Moreover, we calculated the chemical shift of a double-wall tube. We found a diamagnetic shift 
of the isotropic lines corresponding to the atoms of the inner tube due to the effect of the outer 
tube. This shift is in good agreement with recent experiments, and can be easily explained by 
demagnetizing currents circulating the outer tube. 

PACS numbers: 61.46.+w,76.60.Cq,71.15.Mb 



I. INTRODUCTION 

Carbon nanotubes are indubitably among of the most 
interesting materials discovered in the past yeargi^. 
From a theoretical point of view, their quasi one- 
dimensional character imposes quite exotic electronic 
properties. In fact, some time ago it was proposed that 
carbon nanotubes might be good candidates for the elu- 
sive Luttinger liquids^. From a more practical point of 
view, carbon nanotubes turned out to be extraordinary 
materials, as they possess remarkable mechanical and 
electronic properties. This makes them probable build- 
ing blocks for future high-stiffness materials and molecu- 
lar nanodevices^. The electronic properties of single-wall 
carbon nanotubes (SWCNT) can be tuned by varying 
the two integers (n, m) which determine the chirality of 
the tube. For instance nanotubes are either metallic (for 
n = m) or semiconducting. In the latter case, the band 
gap can range from more than 1 eV [e.g., 1.31 eV for the 
(4,3) tube^] to nearly zero [e.g., for tubes with £ — 0, 
where £ = mod {n — m, 3) labels the "family" to which 
the SWCNT belongs]. 

Among various experimental techniques that can be 
used to study and characterize materials, nuclear mag- 
netic resonance (NMR) stands out as being particularly 
suited to study local electronic structureaS'. NMR is 
widely used in structural chemistry and in material sci- 
ence. It measures the local magnetic fields on the nuclei 
generated by the response of the electrons to an exter- 
nal uniform magnetic field. In nonmagnetic insulators 
this response (the chemical shift) is determined by the 
orbital motion of the electrons, and provides information 
about the local atomic environment (coordination and 
geometry) around each nucleus. 

In spite of the usefuUness of the technique, only a few 
groups have up to now used NMR to study carbon nan- 



otubes. The reason is that these experiments are quite 
difficult to perform, due to the presence of magnetic im- 
purities (like Rh, Ft, Ni, etc.) in the samples. These 
impurities, which are added as a catalyst during the 
production of the nanotubes, lead to very broad magic- 
spinning angle (MAS) spectra, with widths of about 30- 
50ppmLiiP. However, we expect this situation to change 
soon, due to new production techniques that use water— 
or alcohol-'^ ^ to enhance the activity and the lifetime of the 
catalysts. In this way, it is possible to produce samples 
with a SWCNT/catalyst weight ratio 100 times higher 
than previous techniques^'', which will clearly allow the 
measurement of NMR spectra to unprecedented accu- 
racy. Another possible approach involves the removal 
of the ferromagnetic impurities in solution. This very 
promising path has been used recently to study the NMR 
response of functionalized (soluble) nanotubesi^. 

In this paper, we perform a ab initio study of the mag- 
netic response properties of infinite carbon nanotubes. 
In particular, we provide benchmark calculations of the 
magnetic susceptibility and NMR ^^C chemical shift of 
semiconducting SWCNT, both isolated and in bundles. 
Furthermore, we present results for double-wall tubes. 



II. COMPUTATIONAL DETAILS 

We use the recently proposed gauge-including projec- 
tor augmented-wave (GIPAW) approach for the calcu- 
lation of NMR response in periodic systemsi2,. This 
approach is based on a linear response formulation of 
density functional theory (DFT), as implemented in the 
plane-wave pseudo-potential code PARATEC^^. Contri- 
butions to the NMR spectra from the core regions are 
taken into account using the projector augmented-wave 
(PAW) electronic structure methodic. We would like to 



2 



remark that the current implementation of the GIPAW 
method only allows the study of semiconducting systems. 
The extension to metals is currently on the way, and first 
results are expected shortly. 

We restrict our study to {n, 0) (zigzag) nanotubes, 
belonging to the two families characterized by £ — 1 
(n = 10, 13, 16, 19) and £ = 2 (n = 8, 11, 14, 17, 20). 
For each of these families, the band gap decreases as l/D, 
where D is the diameter of the tube. (Note that zigzag 
tubes have a relatively small number of atoms in the unit 
cell, which eases the computational burden.) In all cal- 
culations, wave- functions were expanded in plane- waves 
with a cutoff of 40 Ry, and we used the PBE generalized 
gradient approximation to the exchange and correlation 
functional of DFT. For carbon, we used a norm conserv- 
ing TrouUier-Martins^^ pseudo-potential with a core ra- 
dius of 1.75 a. u. The GIPAW reconstruction included 
two projectors per angular momentum for both s and p 
electrons. 

The "isolated" tubes were rolled from a graphene sheet 
with C — C = 1.43 A and were set in a tetragonal super- 
lattice, with a minimum separation of 5.5 A between 
neighboring tubes. For the smaller tubes a sampling of 
20 k-points in the periodic dimension proved sufhcient to 
converge our results to better than 0.25 ppm, while larger 
tubes required 40 k-points at most. We modeled bundles 
by arranging nanotubes on a hexagonal lattice, with an 
inter-tube spacing of 3.3 A. In the plane perpendicular to 
the tube axis, we used a 4 x 4 sampling of the irreducible 
Brillouin zone. Overall, we believe that our numerical 
procedure yields magnetic susceptibilities to a couple of 
percent and isotropic shifts with a precision better than 
1 ppm. 

When calculating the chemical shift tensor, special 
care has to be taken with its macroscopic component 
a{G = 0). This is not a bulk quantity, and its value 
depends on the shape of the macroscopic sample. For 
example, for a spherical or cylindrical sample, and for a 
susceptibility tensor diagonal in the Cartesian axis, this 
quantity can be related to the volumetric macroscopic 
susceptibility xv through the relation (in cgs units) 

(7(0) = -Airaxv , (1) 

where a is a diagonal matrix with values an — 2/3 for a 
sphere and axx = o.yy = 1/2 and azz = 1 for a cylinder 
(in this latter case z denotes the direction of the axis of 
the cylinder). In the case of a bundle of nanotubes, we 
believe that the most correct choice is the cylinder. For 
the isolated tubes the volumetric macroscopic suscepti- 
bility decays as 1/A^, where A is the distance between 
the centers of two adjacent tubes. We are therefore free 
to choose the matrix a in order to speed up the conver- 
gence of the chemical shift with A. We found empirically 
that the choice Uxx — ctyy = 1/2 and azz =2/3 provided 
a faster convergence rate than either the sphere or the 
cylinder. 

Experimental results for the chemical shift are usually 
given relatively to the isotropic shift of tetramethylsi- 



lane (TMS). The reference TMS shielding tensor cttms 
is usually calculated within the same numerical frame- 
work used for the systems under study. This is, e.g., the 
method used in Ref. 17. In this work, we follow a slightly 
different approach. Within our GIPAW method we cal- 
culate the ^'^C shielding tensor of benzene, and use the 
formula 

rTMS _ \so \ I rTMS fr)\ 

"tube — "It^tubc - ^benzene/ "benzene ' \ ) 

where (Ttubo is the calculated shielding tensor of the tube 
and (j'^cnzcnc ^^e average value of the diagonal of the 
calculated shielding tensor of benzene. Finally, <5™2cne 
the experimental isotropic chemical shift in the gas-phase 
of benzene relative to TMS. As the magnetic response of 
the tubes is very similar to the response of benzene, we 
expect in this way to minimize the systematic errors in 
our calculations. The experimental <5j^^f(,jjp was taken to 
be 126.9 ppm.^^, while crjj'enzone calculated at the equilib- 
rium geometry was 40.0 ppm. 



III. MAGNETIC SUSCEPTIBILITY 

The magnetic susceptibility (%) of SWCNTs has 
already been the subject of two theoretical studies: 
(i) Ajiki and AndoiSi using the k.p perturbation method, 
and (ii) Lu2£ using a tigh-binding approach within the 
London approximation. It was found that x increases 
linearly with the nanotube diameter. Furthermore, for 
an applied magnetic field parallel to the tube axis, the 
response X|| was paramagnetic for tubes with £ = 0, while 
for tubes with £ = 1 or 2 the response was diamagnetic. 
In contrast, for a field perpendicular to the axis of the 
tube, all nanotubes were diamagnetic. Finally, it was 
found that the response of the tubes with i = 1 and 
£ — 2 was almost identical. 

The experimental situation is somehow more com- 
plicated, as the experimental samples normally con- 
tain a mixture of tubes of different chiralities, includ- 
ing both metallic and semiconducting tubes. How- 
ever, a recent experimenfeSi estimated the magnitude of 
the magnetic susceptibility anisotropy x\\ ~ Xi- fo be 
~ 14 X 10~^cgs/mol for Inm diameter semiconducting 
tubes, similar to the predicted values of 19 x 10~^ cgs/mol 
(Ref.fT?) and 15 x 10^^ cgs/mol {ReL^. 

In Fig. ^ we plot our ab initio results for the magnetic 
susceptibility tensor of the isolated tubes. In contrast 
with Refs. [T^ and 20 we can clearly distinguish the two 
families of semiconducting zigzag tubes. For the family 
with £ — 2, Xw/D is very weakly dependent on n, while 
for the other family (^ = 1) the absolute value oi x\\/D 
decreases slowly with the diameter of the tube. For x_l 
the situation is reversed: This quantity is almost inde- 
pendent of n for ^ = 1, and its absolute value decreases 
monotonically ior £ = 2. Note that, if we extrapolate our 
results to Z? — > oo, we find values quite consistent with 
Ref. US 
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FIG. 1: (Color online) Magnetic susceptibility of zigzag semi- 
conducting nanotubes, both isolated and in bundles, vs diam- 
eter (D) of the tube. The blue (solid) lines correspond to the 
results of Ref.[li {x\\/D = -5.1 and X±JD = -24.4) and the 
violet (dashed) lines correspond to Ref.l2(]| {x\\/D — —5.4 and 
X±/D = —20.0). The upper axis labels the chirality (n, 0) of 
the tube. All results are in 10~^cgs/(nm mol). 



For a 1 nm diameter tube, we obtain for the magnetic 
susceptibility anisotropy (x|| — X-l) values between 12 
and ~ 18 X lO^^cgs/mol depending on the chirality of 
the tube. This result is in quite good agreement with the 
value x\\ ~ X± = 14 x 10~^cgs/mol measured in Ref. 21. 

Regarding the bundles, we can see from Fig. that 
their magnetic susceptibility is very similar to the one 
of isolated tubes. The main difference is a (fairly rigid) 
decrease of X\\/D by around 0.5 x 10^^ cgs/ (nmmol) and 
of X±/D by around 1 x 10~® cgs/ (nmmol). 



IV. CHEMICAL SHIFT 

Only a few groups have studied carbon nanotubes via 
NMR, whether experimentally or theoretically. Experi- 
mental ^^C MAS spectra of pristine tubes yield isotropic 
shiftsSI (5iso of 126ppm'''^, 124 ppm.^^, and 116 ppm^. The 
full width at half maximum (FWHM) of the isotropic 
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FIG. 2: (Color online) Isotropic chemical shifts for the iso- 
lated zigzag nanotubes, versus diameter (D) of the tube. The 
hnes are the fits to the numerical results given by Eq. 



peak is fairly large, being around 50ppm'''® or BOppmSt. 
This width was interpreted as a distribution of isotropic 
lines, resulting in part from the presence of both semi- 
conducting and metallic tubes in the sample. The static 
spectra are quite broad, spawning more than 200 ppm. 
It was also noted that the presence of residual Ni and Co 
magnetic particles can contribute to the further broad- 
ening of the spectra. 

The principal values of the chemical shift tensor can be 
obtained from the NMR spectra through a fitting pro- 
cedure. In this way, Ref. found S™^ = 195 ppm, 
Sji^^ = 160 ppm, and 6ji^^ = 17 ppm (which leads to 
A(5anis — — 160.5ppm and 77c = 0.33), while in Ref. 
it was obtained S™^ = 240 ppm, S™^ = 171 ppm, and 
5™S ^ _36ppni (A<5anis = -241.5ppm and r/c = 0.43). 
The large difference between the two experiments (see 
Fig. 0J can be understood from the difhculty of obtain- 
ing good samples. 

On the theoretical side, the first attempt at describ- 
ing the NMR response of carbon nanotubes^^ was per- 
formed within a tight-binding scheme and using the Lon- 
don approximation (i.e., retaining only the London ring- 
current contribution to the chemical shift and neglecting 
the Pople corrections^). To justify this approximation 
it was argued that the Pople contribution gave rise to 
a global shift, independent of chirality, that could not 
be used to distinguish between different tubes. Within 
this framework, it was estimated that the difference in 
the isotopic chemical shift between metallic and semi- 
conducting tubes amounted to 12 ppm. (Absolute shifts 
were not given due to the neglection of the Pople cor- 
rection.) Furthermore, it was found that the NMR of 
both metallic and semiconducting is mostly insensitive 
to diameter and chirality. 

More recently the ^'^C chemical shift of finite-size (9,0) 
single wall carbon nanotubes was investigatedii. The 
tube fragments studied were capped by either C30 or hy- 
drogen. For the longest fragment it was found <5™^ ~ 
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FIG. 3: (Color online) Chemical shifts relative to TMS for the zigzag nanotubes in bundles versus radius of the tube. The blue 
(solid) lines connect the average values for each tube. The green (dashed) line is the results for the isolated tubes presented 
before. 



133 ppm for carbon atoms at the center of tube, in rea- 
sonably good agreement with experimental data. For the 
CsQ-capped tubes, apparently the best models for the in- 
finite systems, the shifts in the middle of the fragment 
were invariably smaller than those at the end (where they 
reached S™^ « 155 ppm). 

A. Isotropic shift 

Our results for the isolated tubes are presented in 
Fig. 12 The two families can be clearly identified: The 
family with £ = 1 has a slightly larger chemical shift than 
the family with i = 2. The results can be fitted by the 
relation 

cTMS _ / 10.70/2? + 116.0 if ^ = 1 

\ 9.83/L»+ 116.0 if ^ = 2 ' ^ ' 

where D is the diameter of the nanotube in nanometers, 
and (5™^ is the isotropic chemical shift relative to TMS 
in ppm. According to this formula, in the limit of infinite 
radius the isotropic chemical shift tends to 116.0 ppm. 
This value compares reasonably well with the estimate 
of 128 ppm for ^Jraphonc '^^ * perfect graphene layer (with 
zero density of states at the Fermi level)^*^^. For a nan- 
otube with a diameter of 1.2 nm, like the ones common 
in experimental samples, we obtain a Siso ~ 125 ppm, in 
excellent agreement with experimenliiSiSiS^. Note that, 
as already indicated in Refs.0 and|^ the NMR chemical 
shift depends very weakly on the diameter and on the 
chirality of the semiconducting nanotubes: the isotropic 



shifts of the zigzag tubes studied in this work lie between 
122 ppm for the smallest tube [the (8,0)] to 132 ppm for 
the (20, 0) tube. 

The isotropic shifts of the bundles are represented by 
red crosses in Fig. |31 Since the hexagonal lattice breaks 
the chiral symmetry of the tubes, the carbon atoms 
around the circumference of the tube are no longer equiv- 
alent, resulting in a slight spread of S'^^^ in the unit 
cell. For the nanotubes studied here, this spread of the 
isotropic shift was quite small (around 3 ppm), and fairly 
independent of the radius of the tube. It can also be seen 
that the average shift in a bundle is nearly equal to the 
isotropic value of the isolated tube^. 

We are now in condition to discuss the large width 
MAS spectra, that amounts to 30-50 ppm. Several effects 
can contribute to this width: 

(i) Experimental samples are composed of nanotubes 
of different diameters and chiralities. Assuming that the 
spread of diameters is around 0.5 nm, and the distribu- 
tion is centered around 1.2 nm, from Fig [21 we obtain 
a spread of around 4 ppm. Furthermore, the difference 
in chemical shift between semiconducting and metallic 
tubes was estimate to be 12ppm&. 

(ii) As experiments are usually performed with samples 
of SWCNTs arranged in bundles, we have to take into 
account the spread of isotropic chemical shifts due to the 
breaking of chiral symmetry. As we saw from Fig. lathis 
amounts to around 3 ppm for zigzag tubes. 

(iii) As mentioned before, in order to obtain the macro- 
scopic component of the chemical shift tensor, we as- 
sumed that the bundle was of cylindrical shape. However, 
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FIG. 4: (Color online) Experimental versus simulated static 
spectra for the carbon nanotubes. The curve labeled 'exp 1' 
was taken from Ref. .22s while 'exp 2' was taken from Ref. 
The theoretical curve corresponds to the static spectrum of a 
bundle of (14,0) nanotubes, convoluted with a Lorentzian of 
30 ppm FWHM. 



in a real sample, the bundle can assume shapes quite dif- 
ferent from a perfect cylinder. To address this issue, we 
performed calculations assuming a spherical shape for the 
bundle. For the (14,0) tube, a typical nanotube of 1.1 nm 
diameter, the isotropic shift calculated for the cylinder 
was 125.3 ppm, while the sphere gave 122.4 ppm. Note 
that, even if the difference, 3 ppm, is much smaller than 
the experimental width of the MAS spectrum, it has nev- 
ertheless the same magnitude as the difference between 
nanotubes of similar diameter. 

(iv) Finite-size SWCNTs can be either open, or capped 
by fuUerene-like structures. It was found in Ref. J7 that 
the carbon atoms close to the ends of these (finite) tubes 
can have an isotropic shift up to 20 ppm larger than the 
atoms far from the ends. However, although not infinite, 
the SWCNTs present in experiment are quite long, so we 
do not expect this effect to contribute significantly to the 
NMR spectra. 

(v) Other effects: Clearly every experimental sample of 
SWCNTs contains defects and magnetic impurities which 
will contribute to the broadening of the MAS peak. Fur- 
thermore, there may be interactions of the semiconduct- 
ing nanotubes with the conduction electrons of metallic 
tubes contained in the bundle, etc. 

We can see that effects (i)-(iii) can contribute around 
^ 20 ppm to the width of the MAS peak. We therefore 
expect that the remaining ~10~30ppm of the experimen- 
tal widt h'^1^1^ are due to (v). This means that by using 
better samples with a much lower concentration of mag- 
netic impuritiesifliii the resolution of ^"^C NMR spectra 
can be improved by at most a factor of 2. 

B. Anisotropic shifts 

Also in Fig. |21 we show the principal values of the 
chemical shift tensor for the bundles. The first princi- 



pal value, 5™^ corresponds to the longitudinal (along 
the tube axis) direction, (5™^ lies on the orthoradial di- 
rection, and is the radial component of the chem- 
ical shift tensor. As expected, the chemical shift is ex- 
tremely anisotropic, with the longitudinal and orthora- 
dial directions corresponding to the in-plane directions 
of a graphene sheet (or of a benzene molecule), and the 
radial component corresponding to the direction perpen- 
dicular to the sheet. As for the isotropic shift, the in- 
dividual principal values are spread due to the interac- 
tion between the tubes in the bundle, but this time by 
as much as 10 ppm. Note that Ref. found a value of 
20 ppm for the broadening of the tensor due to the in- 
fluence of the first neighbor tube (and quasi-independent 
on the electronic properties of this first neighbor) . There 
is, therefore, a partial cancellation when we include all 
the six neighbors in the hexagonal lattice. Note that this 
10 ppm spread is much larger than the 3 ppm spread of 
the isotropic line. It is also much larger than the dif- 
ference between the two families of nanotubes we have 
studied. 

In Fig. 0] we compare a simulation of the static spec- 
trum of a (14, 0) tube with experimental results. The the- 
oretical spectrum shows a typical powder line-shape with 
principal values 5™^ = 198 ppm, (5™^ = 188 ppm, and 
^ -iippm (A(5a„is = -204 ppm and r;c = 0.07). 
The two experimental curves are quite different from each 
other, so care must be taken when comparing theoretical 
to experimental spectra. In general, we believe that our 
results are in quite good agreement with experiments. 
The main difference turns out to be the value for the 
asymmetry factor, which is considerably smaller in our 
calculations than in the experimental works. 



V. DOUBLE- WALL TUBES 

Recently, Simon et al. studied double-wall nanotubes 
using NMR22,. Their samples were produced through the 
high- temperature annealing of isotope enriched fuUerenes 
encapsulated in single-wall tubes. In this way, they were 
able to measure directly the chemical shift of the inner 
tube. They found an isotropic shift of 111 ppm, which 
is significantly smaller than the isotropic shift of larger 
single wall tubes. The authors pointed out two possible 
reasons for this difference: (i) the stronger curvature of 
the inner tube, or (ii) diamagnetic shielding due to the 
outer tube. Note that in our calculations the chemical 
shift increases with decreasing radius (see Fig. ^ , which 
rules out (i). 

To study this problem, we performed ab initio cal- 
culations of the isolated double-wall tube (8,0)@(16,0). 
We obtained for the isotropic chemical shift of the atoms 
composing the inner tube a value of 111 ppm [to be com- 
pared with 132 ppm for the isolated (8,0) tube], implying 
a diamagnetic shift due to the outer tube of 21 ppm. On 
the other hand, the outer tube had an isotropic chemical 
shift of 126 ppm [to be compared to 124 ppm of the iso- 



6 



lated (16,0) tube]. Note that our ab initio result for the 
inner tube (lllppm) coincides with the value measured 
in Ref.S 

In order to understand our results, and to study the 
effect of diamagnetic shielding, we will resort to a simple 
classical model. Let us imagine that the outer nanotube 
has the form of a long hollow cylinder with axis parallel 
to z. We will distinguish two cases, when the probing 
magnetic field Bq \\ z and when Bq _L i. In the former 
case, the induced current j on the surface of the cylinder 
can be written (in linear response) as j = (3\\BaO-, where 
/3|| is a constant. The magnetization per unit length can 
then be calculated from (in cgs units) 



M = ^— (Sr r X iir) = 
2cL J "'^ ' 



(4) 



where c denotes the speed of light, L is the lattice pa- 
rameter in the periodic direction, D is the diameter of 
the cylinder, and z is the unit vector of direction z. The 
volume of integration is one unit cell. Remembering that 
M = j^ncX\\Bo, where X|| is the molar susceptibility 
parallel to the axis of the cylinder, A^a is the Avogadro 
number, and nc is the density of carbon atoms per unit 
area of the nanotube (nc — 37.71 nm^^), we can obtain 
the induced current 



(5) 



On the other hand, it is a simple excercise of magneto- 
statics to calculate the magnetic field generated by Bq 
and by the induced current. The total magnetic field 
(Bin) inside the cylinder (i.e., the field felt by the inner 
nanotube) is constant. Furthermore 



Bin - Bo 4nc 

= 47r^^— XII 



B, 



NaD' 



(6) 



Note that as xy|| is negative for our nanotubes, Eq. © 
gives rise to a diamagnetic shift. The same derivation can 
be done for an external field perpendicular to the tube 
axis. In this case, the current is, for a magnetization in 
the y direction, j = cos{9)z on the sides of the cylinder, 
and j = ±(3±Box on its bases. The total (external plus 
induced) field inside the tube is still constant, and it reads 



Bin — Bq 

Bo 



(7) 



We can now calculate the diamagnetic shift in the inner 
tube due to the presence of the outer tube. As a typi- 
cal example we chose the double-wall tube (8,0)@(16,0). 
Inserting its diameter (1.26 nm) and the calculated sus- 
ceptibility (x|| = —9.60 X 10~^cgs/mol and x± — 
-24.86 X IQ-^cgs/mol) of the (16,0) tube in Eqs. ® 
and O we obtain that the shift due to the shielding is 
^shield ^ _24ppm and Sff'"^'^ = S^^'"^'^ = -31ppm. This 
leads to an isotropic shift of 5fs'j,'°''^ = — 29ppm. From our 
calculations, the isotropic shift of the isolated (8,0) tube 
is 131 ppm, which leads to a total chemical shift for the 
inner tube of ~ 102 ppm. This is a fairly good agreement 
with our ab initio results for such a crude model, and 
fully justifies its use in these systems. 



VI. CONCLUSIONS 

We performed calculations of the magnetic suscepti- 
bility and ^'^C chemical shift for semiconducting zigzag 
carbon nanotubes, both isolated and in bundles. Our re- 
sults are in very good agreement with experiment. For 
all nanotubes studied, the isotropic shift lies between 
120 ppm and 135 ppm. This interval of 15 ppm is quite 
narrow, and is at the limit of the precision that NMR 
experiments have obtained for these systems. However, 
we expect that, with the advent of better samples, NMR 
will become a key tool to study and characterize carbon 
nanotubes. 
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